Frozen into stripes: fate of the critical Ising model after a quench 
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In this work we study numerically the final state of the two dimensional ferromagnetic critical Ising 
model after a quench to zero temperature. Beginning from equilibrium at Tc, the system can be 
blocked in a variety of infinitely long lived stripe states in addition to the ground state. Similar 
results have already been obtained for an infinite temperature initial condition and an interesting 
connection to exact percolation crossing probabilities has emerged. Here we complete this picture 
by providing a new example of stripe states precisely related to initial crossing probabilities for 
various boundary conditions. We thus show that this is not specific to percolation but rather that 
it depends on the properties of spanning clusters in the initial state. 
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Introduction — Somewhat surprisingly, even with such 
a simple set up as the two dimensional nearest neighbour 
ferromagnetic Ising model quenched to zero temperature, 
a variety of final states can be reached with non trivial 
probabilities. It was first noticed by Lipowski [T] that 
an anomalous scaling for the equilibration time of the 
kinetic Ising model arises from the existence of stripe 
states. In an interesting series of papers, Krapivsky, 
Redner and collaborators noticed as well the presence of 
stripes states after a quench from infinite temperature, 
and analysed both the dynamics and the final states. 
They measured numerically the probability of getting 
stuck in stripe states and found that it was around 1/3 
but were at first unable to explain it [21 [3] . In the same 
situation but with different goals, Sicilia et al. studied 
the geometry of spin clusters, and noticed that after a 
few Monte-Carlo steps at T < Tc, an infinite temper- 
ature system on a square lattice is very similar to the 
critical percolation point for the spin clusters jUIS]. So, 
even though the initial occupation probability of 1/2 for 
a type of spins is inferior to the critical site percolation 
probability for the square lattice p^ ~ 0.59, the system is 
rapidly very similar to a critical percolation system. This 
correspondence enabled Krapivsky and collaborators to 
give an explanation for the probabilities of appearance 
of stripe states in terms of critical percolation probabili- 
ties, be it for free or periodic boundary conditions [SlII]. 
The study of these quantities has a long history start- 
ing with Cardy and his eponymous formula 0. This 
formula gives the probability of the existence of an in- 
cipient spanning cluster at the critical percolation point 
in terms of hypergeometric functions. Since then a num- 
ber of percolation probabilities have been studied. One 
aspect that is particularly interesting is that although 
Cardy and others found those formulas using conformal 
field theory (CFT), i.e. in a non rigorous way mathemat- 
ically speaking, mathematicians have developed rigorous 
tools to tackle those systems. Schramm introduced the 
stochastic Loewner evolution (SLE) [5] which describes 



numerous physically occurring curves, and this lead to 
many results (see jlOj with Lawler and Werner e.g.). 
Around the same time, Smirnov proved rigorously Cardy 
formula and the conformality of the critical percolation 
point for site percolation on the triangular lattice [IT] . 

The correspondence discovered by Krapivsky et al. is 
thus very interesting since it relates a non-equilibrium 
situation to exact results. It is remarkably accurate and 
based on the observation of Sicilia et al. and general con- 
siderations on coarsening dynamics and it was supported 
by numerical results on a lattice with different aspect ra- 
tios. The underlying question to those studies is to what 
extent the initial condition infiuences the final one? In 
this regard, it is interesting to study the same situation 
with a different initial condition to extract the universal 
behaviors from the rest. For the Ising model, the results 
discussed above have to hold for any starting tempera- 
ture above the critical temperature since for an initial 
condition at T > Tc the long distance properties are gov- 
erned by the infinite temperature fixed point. Since the 
subcritical region is trivial from the spin clusters point of 
view, the critical temperature point only remains and is 
absolutely non trivial concerning spin clusters as we will 
disscuss below. Actually, the persistence of the initial 
condition is also very interesting to study before reach- 
ing the final state, i.e. during the equilibration. Several 
works have dealt with the issue of cluster dynamics after 
a quench either to T < Tc [H [S] or to Tc [H] for the Ising 
model and for the Potts model [131 [H] and were able to 
clarify the influence of the initial condition. Indeed those 
works show that the initial properties of spin clusters are 
retained at distances bigger than a dynamical length- 
scale ^{t) ~ t^'^ using dynamical scaling arguments and 
numerical checks. 

In the following we study the ferromagnetic Ising 
model whose hamiltonian is written as: 

'H = -Y,JS.^S,, (1) 

where the sum is over all nearest neighbors pairs of spins 




FIG. 1: (color online) Crossing probabilities for free bound- 
ary conditions versus r. The lines corresponding to 7r^,7rJ 
and ttJjj, are from the results given in [T^ for the critical Ising 
model. The points correspond to our data for tt^, ttJ) and 7r}J„ 
obtained after a quench at zero temperature. 



of a two dimensional system and J > 0. The system 
considered is rectangular and three boundary conditions 
will be considered : free, periodic and fixed. The lattices 
used will be discussed in each case. 

Method — In the simulation that we will consider here, 
we started from a system equilibrated at the critical tem- 
perature. The equilibration is obtained with the Wolff 
cluster algorithm [IB^. We generate a large number of 
independent configurations, at least one million for each 
size considered. For each configuration, we determine the 
probability that there exists at least one cluster of spins 
crossing the horizontal or vertical direction. We denote 
these probabilities 7r£ or tt^. The probability that a spin 
cluster is crossing in both directions will be denoted 7r£^ . 
Next, we perform a quench at T = 0. We then let the 
system evolve with a Glauber dynamics with a suitable 
algorithm described below until it reaches a stable state. 
We record similar (stripes) probabilities that we denote 
by TT^ , ttJJ for the horizontal or vertical stripes and ttJJ^ for 
the stripes in both directions, i.e. a ground state. 

One difficulty is that it can take a tremendous amount 
of (Monte-Carlo) time to reach the final state where no 
spins can be updated without increasing the energy of 
the system because the system tends to wander in long- 
lived iso-energetic states. A simple kinetic Monte Carlo 
algorithm [TTj with Glauber dynamics bypasses this issue 
quite easily. It essentially consists in an efhcient sampling 
of the updateable spins which greatly reduces computa- 
tion time to reach a final state. In a T = dynamics, 
the only possible transition rates on even-coordinated lat- 
tices are 0, 0.5 and 1 respectively for spins with strictly 
more than half of their neighbours of the same colour as 
theirs, exactly half and strictly less than half. Let us call 



n+, no and n_ the number of spins in this three cate- 
gories and R the sum of all possible rates, with obviously 
R = n^ + no/2. A list of the n_ + rig updateable sites is 
created. Now, to update the system, a random number 
u s]0, 1] is drawn, and if uR < ri-, one of the ri_ spins is 
flipped and otherwise one among the uq ones. To update 
the time i, we draw another random number u' s]0, 1] 
and t -i— t-|-ln(l/u')/it!. Then the list of updateable sites 
is regenerated and the system can be updated again until 
n^ + Hf) = 0. 

Free Boundary conditions — We will first describe our 
results for a system with free boundary conditions. The 
probabilities tt^, tt^ and 7r£^ were already considered 
in equilibrium at T^ by Lapalme and Saint- Aubin [15] . 
These authors measured these quantities on the triangu- 
lar lattice. They also tried to find a way to predict the 
behaviour of these crossings in a way analogous to the 
Cardy formula for percolation. They obtained a differen- 
tial equation on a four point correlation function related 
to 7r£ but were unable to solve it so they resorted to a 
numerical solution in good agreement with their mea- 
surements. 

We considered systems of size L x L/r with L — 256 
while varying r. We only study the case r < 1, the case 
with r > 1 being obtained by a trivial duality relation 
exchanging the vertical and horizontal directions. Our 
results are shown in Fig. [Tl In this figure we compare 
the situation at Tc and the final states of the evolution 
at T = 0. To do so we present the probabilities of get- 
ting a crossing 7r°, 7r° and 7r°^ for various values of ratio 
r obtained at T = as crosses and in the critical case, 
tt^jTtJJ and 7r^„ shown as lines. The agreement between 
the tt'^ and the tt" is excellent, this is one of the main re- 
sults of the present work. The values of the tt'^ are those 
obtained by Lapalme and Saint- Aubin in [15| . Note that 
these authors considered the case of a triangular lattice 
while our measurements have been done on a square lat- 
tice. As a check, we also repeated the measurement of 
the TT^ on the square lattice and of the tt" on the tri- 
angular lattice and our results are in perfect agreement 
with the ones of Lapalme and Saint- Aubin obtained on 
the triangular lattice. This supports the universality of 
these quantities. 

In conclusion, the crossing probabilities are the same 
for the equilibrium at Tc and in the blocked states ob- 
tained after a quench at T = starting from the critical 
point. This first result confirms the idea that the final 
state is dictated by the topology of spin clusters in the 
initial one. 

Periodic Boundary conditions — To make this fact 
more explicit we move on to the case of periodic boundary 
condition (PBC). Here the situation is a bit more com- 
plicated since clusters can wind in various ways around 
the torus. We can nonetheless define as previously the 
percolation probabilities 7r£, tt^, tt^^ and 7r°, 7r°, 7r°^ of 
interest to us. We consider the triangular lattice. With 




FIG. 2: (color online) Two possible types of stripes of the al- 
ternate triangular lattice with periodic boundary conditions. 
The vertical stripes on the left can exist for arbitrary aspect 
ratios r contrarily to the stripes on the right which only exists 
for r — 2n/\/3 with n £ N*. Actually the stripe on the right 
wraps in both directions around the torus. 



this choice of lattice and boundary condition only verti- 
cal stripes are stable at T = for arbitrary aspects ratios 
(see Fig. [2|. This means that tt^ = and tt^^ = 1 ^ ^iS 
thus we only need to consider ttJJ. As is shown in Fig. pi 
diagonal stripes are stable in addition to the vertical ones 
for the aspect ratios r = 2n/y/2> with n S N*, but their 
rarity (TTdiag is at most 6.10"*) forbids their study nu- 
merically as has been done in [7]. The factor -s/S/S in the 
definition of the aspect ratio takes into account the actual 
horizontal system size. We have seen that different perco- 
lation probabilities of same spin cluster for free boundary 
conditions has been studied numerically and analytically 
in [15 . We found no such study for the PBC. This mo- 
tivated us to extend the work of Arguin [TH] which deals 
with the probability that a given configuration contains 
a Fortuin-Kasteleyn cluster winding a times horizontally 
and b times vertically around the torus. We have been 
able to obtain an explicit formula for such a probabil- 
ity for an Ising spin cluster [19] . In the case where the 
spin cluster winds only vertically around the torus, the 
formula reduces to: 



<{r) = 



1 e,{^)-e,{%)-29,{'^) 



2r|77(^r)||02(^r)| 



(zr) I + |04 [ir) 



(2) 



where rj is the Dedekind 77 function and the 9i are the Ja- 
cobi 9 functions. In Fig. [3} the results for the triangular 
lattice with PBC are presented. The agreement between 
TT^ , TT^ (dots) obtained from the simulations and the the- 
oretical prediction from Eq. ^ for tt^ (dashed line) is 
very good. This second case confirms the link between 
probabilities of blocked stripe states and initial crossing 
probabilities of spin clusters. 

Fixed Boundary conditions — Finally, we will look at 
a last case in which wc impose fixed boundary conditions 
with S" = — 1 on the left and right sides of a rectangle 
and with 5 = +1 on the top and bottom sides. This 
case is interesting because it has been extensively stud- 
ied analytically at Tc by several authors using CFT and 




FIG. 3: (color online) Vertical percolation probability tt^ at 
Tc (red crosses) and vertical stripe probability ttJ) at T = 
(blue squares) for various aspect ratio r = Ly/{\/iLx/2) with 
Ly — 128. The dashed line corresponds to the expression in 
Eq. ||2|. 



multiple stochastic Loewner evolutions [20H21]) so there 
exists a number of theoretical predictions to compare our 
simulations to; moreover it is really easy to simulate and 
analyze as we will see below. 

These boundary conditions will force the existence of 
two interfaces. The two possible configuration types are 
represented in Fig. [2] In the first case, the first interface 





FIG. 4: Schematic representation of the two configuration 
types with fixed boundary conditions. Is is easy to see that 
beginning with the left (respectively right) configuration we 
will end up in a state where all spins equal +1 (—1) except 
the fixed ones. 

goes from the left top corner to the left bottom corner 
and the second interface connects the right top corner to 
the bottom right corner. In the second case the inter- 
face goes from the left top corner to the right top corner 
and the second interface connects the bottom left corner 
to the bottom right corner. The probabilities for these 
two situations at T^ can be written in several equivalent 
forms, we use Kozdron's version written only in terms of 
the hypergeometric function F |22| . The probability for 
the left situation in Fig.|4]is given by: 



P{x) 



^ V -J J "J; 3 j 1 ^) 



-f* ( T : •Jj -J j 2; j + -T ( „ , O, „ ; 1 X) 



(3) 



where x € (0, 1) is related to the aspect ratio r by: 

K{x) 
^ " K{l-x) 



(4) 



with K{x) the complete elliptic integral of the first kind. 
For the case of infinite initial temperature, the probabil- 
ities for the two situations of Fig. |4] are given by Cardy 
formula [5j. Indeed, in the case of percolation, the fixed 
boundary condition considered in [5DH12] coincides with 
the situation considered by Cardy. 

As previously we quench the system from Tc to T — Q 
but also from T = oo to T = as we found no mention 
of this case in the literature. With this boundary con- 
dition the analysis of the final state is really easy, the 
sign of magnetization suffices to indicate the state of the 
system as bulk spins are all of the same sign in the end. 
We show the results of our simulations for L — 1280 for 
the initial temperatures T^ and T = oo in Fig. [5] The 
agreement with the theoretical prediction (dashed lines) 
is again excellent in both cases. 
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FIG. 5: (color online) Crossing probabilities versus r. The 
dashed lines correspond to the formula proved in |20H22) for 
the critical Ising model and in 8 for percolation. The symbols 
correspond to our data for vr^' (blue filled square) (respectively 
■K^'^ (green dots)) obtained after a quench at zero temperature 
from Tc (respectively T — oo). Note that the values measured 
for ttJ (red crosses) are indistinguishable from the probability 
n° since they are numerically very close. 

Conclusions — We have extended the connection be- 
tween crossing probabilities and the probabilities of exis- 
tence of stripe states after a T = Glauber dynamics to 
the case of the critical Ising model for various boundary 
conditions. We have obtained clear results showing that 
the final state of the evolution is strongly correlated to 



the initial condition, similarly to what was already ob- 
served in [21 121 m [7] for the quench from T = oo toT = 
in a similar setup. We expect that this is a general fea- 
ture and it would be interesting to check the generali- 
sation of this result for other models. In a preliminary 
check, we have obtained similar results for the Q — 3 
Potts model f^. We leave for the future similar studies 
for more complicated models or even the case of larger 
dimensions. 
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